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ABSTRACT 

We explain the Ha emission of the cometary knots in the Helix Nebula (NGC 7293) with an analytical model 
that describes the emission of the head of the globules as a photoevaporated flow produced by the incident ionizing 
radiation of the central star. We compare these models with the Ha emission obtained from the HST (Hubble 
Space Telescope) archival images of the Helix Nebula. From a comparison of the Ha emission with the predictions 
of the analytical model we obtain a rate of ionizing photons from the central star of about 5 x 10 45 s" 1 , which is 
consistent with estimates based on the total H/3 flux of the nebula. We also model the tails of the cometary knots 
as a photoevaporated wind from a neutral shadow region produced by the diffuse ionizing photon field of the 
nebula. A comparison with the HST images allows us to obtain a direct determination of the value of the diffuse 
ionizing flux. We compare the ratio of diffuse to direct stellar flux as a function of radius inside an HII region 
with those obtained from the observational data through the analytical tail and head wind model. The agreement 
of this model with the values determined from the observations of the knots is excellent. 
Subject headings: ISM: structure — planetary nebulae: individual (NGC 7293) — stars: AGB and post-AGB 



1. INTRODUCTION 

The small scale structure of the Planetary Nebula (PN) 
'. known as the Helix Nebula (NGC 7293, PK 36-57° 1) is charac- 
terized by many thousands of small knots. These knots (which 
were first reported by Vorontsov-Velyaminov 1968) have a 
cometary shape, with their tails pointing away from the cen- 
tral source. Groundbased work (Meaburn et al. 1992, 1996, 
1998) and Wide Field and Planetary Camera HST observations 
(O'Dell & Handron 1996; O'Dell & Burkert 1997; Burkert & 
O'Dell 1998) revealed a multitude of spatially resolved knots, 
about 3500, with highly symmetric appearance. The knots have 
also been detected in CO (Huggins et al. 1992) and, at least out- 
side the main nebula, in C I (Young et al. 1997). 

This PN is one of the closest to us with a parallax distance of 
213 pc (Harris et al. 1997). Other methods have been applied to 
determine the distance to this PN, giving values ranging from 
120 to 400 pc. Optical images show a complicated morphol- 
ogy characterized by a helical structure in Ha and [N II] and a 
more elliptical shape in [O III]. The deprojection of these im- 
ages has proven to be difficult, with suggestions of both ellipti- 
cal/toroidal shapes and bipolar ones (Meaburn & White 1982) 
being made. O'Dell (1998) has suggested that the ionized re- 
gion is a disk. 

The central star is very hot and has a low luminosity, indi- 
cating that it is well down the cooling track of its post-AGB 
evolution. The temperature of the central star of the Helix has 
been measured by the H/3 Zanstra method to be ~ 1.2 x 10 5 K 
(Gorny, Stasiriska & Tylenda 1997). 

The high effective temperature combined with the fairly low 
luminosity and the large size of the PN , about 0.5 pc in diam- 
eter at a distance of 213 pc, leads to a complicated ionization 
structure with co-existing regions of high and low ionization 
character (O'Dell 1998; Henry, Kwitter & Dufour 1999). The 
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low luminosity of the star also implies a very low density wind, 
making it virtually undetectable. However, with its expected 
high velocity (Vqo ~ 6000 km s" 1 if radiation driven) it could 
still have a dynamical effect on the PN structure. 

The measured densities in the nebula are low, about 60 cm" 3 
(O'Dell 1998), in contrast with the densities derived for the 
neutral/molecular globules, n^oi ~ 10 6 cm" 3 (O'Dell & Han- 
dron 1996). This leads to the curious result that a substantial 
fraction of the mass of the PN is in the form of these neutral 
globules. The aspherical shape of the PN is also found in the 
distribution of the knots and in the large scale CO emission. 
Meaburn et al. (1998) have shown that the knots follow the 
same velocity distribution as the CO, although at a lower typical 
expansion velocity. Young et al. (1999) found the deprojected 
expansion velocities of the knots and the CO ring to be 19 and 
29 km s" 1 , respectively. 

Some other authors have treated the photoevaporation of 
clumps in an ionizing radiation field. Bertoldi & McKee (1990) 
developed an approximate analytic theory of the evolution of a 
photoevaporating cloud exposed to the ionizing radiation of a 
newly formed star, finding an equilibrium cometary cloud con- 
figuration. Johnstone, Hollenbach & Bally (1998) modelled 
the photoevaporation of dense clumps of gas by an external 
source of ultraviolet radiation including thermal and dynami- 
cal effects. Mellema et al. (1998) study the evolution of dense 
neutral clumps located in the outer parts of a planetary nebula 
(like the Helix nebula). 

2. ANALYTICAL MODEL FOR THE PHOTOEVAPORATED FLUX OF 
THE COMETARY KNOTS 

Let us discuss a simple, analytical model for cometary knots 
embedded in a Planetary Nebula (PN). First, we derive some 
expressions that relate the ionizing stellar flux and the particle 
flux of the photoevaporated wind in a similar way to Mellema 
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Photoevaporating flows from the cometary knots in the Helix nebula (NGC 7293) 




FIG. 1 . — Schematic diagram of the cometary head. F/, is the incident flux at the surface of the neutral knot, and R/, and are the radius and the density of the 
neutral knot (respectively). 



et al. (1998) to find the Ha emission for the cometary heads. 
Canto et al. (1998) proposed the idea of the tails as neutral 
shadow regions behind the clumps, and studied the complex 
time-evolution of the resulting flow. In this paper, we model the 
tails behind the Helix clumps as a cylinder of neutral material 
being photoionized by the diffuse ionizing flux of the nebula. 

2.1. A model for the cometary head 

Let us consider the problem of a hemispherical, neutral knot 
of radius R/, which is being photoionized by the radiative field 
emitted by the central star of the ionized nebula. A schematic 
diagram of this configuration is shown in Figure [I]. 

A fraction of the stellar flux arrives at the knot surface, pho- 
toionizing the neutral knot material and feeding a photoevapo- 
rated flow. The remaining ionizing photons are absorbed in this 
photoevaporated flow. The neutral knot is being photoevapo- 
rated, but we consider a quasi-steady state in which the radius 
of the knot R/, is almost constant with time. The incident flux at 
the knot surface is : 

F h (6) = n h (9)v(R h ), (1) 

where F/, is the flux of ionizing photons incident normal to the 
knot surface (per unit area and time), «/, and v(i?&) are the den- 
sity and the velocity of the photoevaporated wind at the knot 
surface . We note that the rate of photoionizations depends 
on the angle of incidence 9 (see Figure [j]), having a maximum 
value in the front of the cometary head. 

In order to obtain the density profile, we consider the conser- 
vation of particles in the photoevaporated wind : 

R 2 n(R,0MR) = Rln h (8)v(R h ), (2) 

where R is the radius directed outwards from the knot. 
We assume that the ionization front is thin compared with the 
effective thickness h of the ionized flow and that h in turn is 
small compared with /?/,. In this case, if all diffuse photons are 
absorbed "on-the-spot", the photoionization balance is given by 



F* cos(#) ~ F h (9) + h n\{9) a B , (3) 

where as is the case B recombination coefficient of hydrogen, 
which is assumed to be constant with position. The flow thick- 
ness h is defined by : 

h = ujR h , (4) 
where the parameter ui is defined through the relation : 

POO 

ujn 2 h R h = I n 2 (r)dr. (5) 
JRh 

This parameter ui depends on the density profile of the photo- 
evaporated wind. For a D-critical ionization front (v(/?/,) = c,, 
where c, is the isothermal sound speed of the ionized gas), 
uj ~ 0.1 (Henney & Arthur 1998). From equations (1-5) we 
then obtain a relation between the ionizing flux from the cen- 
tral star F*, and the incident flux at the knot surface Fi,(9) : 

^cos(9)^F h (9) + F 2 (9)^^ 1 (6) 

The term on the left is the ionizing stellar flux as a function of 
the angle of incidence, the first term on the right is the incident 
ionizing flux on the knot surface that produces the photoion- 
izations of the neutral material and the second term is the flux 
absorbed in the photoevaporated wind. If we define : 



Rh u OL B 

equation (6) can be written as : 

F, CQ ^9)-F h {9) + ^P-, (8) 

We can see from equation (8) that it is possible to find the 
incident flux F/,{9) as a function of the ionizing flux solving 
the quadratic equation to obtain : 
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(9) 



We find two different limiting situations : 

If £/,/ cos(0) «F* we have the "Recombination Dominated 
Regime " in which the fluxes are related by the expression : 



F h (9) « v^ivcos^), 



(10) 



In this case, an important fraction of the ionizing flux is ab- 
sorbed in the photoevaporated wind and only a small fraction 
of the stellar flux arrives at the knot surface. The recombina- 
tions in the column beyond Rf, mainly balance the stellar flux. 

If £/,/ cos(6*) >>F t we have the "Flux Dominated Regime ", 
in which : 



F h (8)~F+ cos(0), 



(11) 



In this regime almost all the stellar flux arrives at the knot sur- 
face and there is no absorption in the photoevaporated flow. The 
paticle flux at R/, is roughly equal to the stellar flux. 

In order to compute the total rate of Ha photons emitted by 
the head of a cometary knot, assuming wClwe have to eval- 
uate the integral : 



Sho 



■k/2 



2 7T R\ qho n\(6) sin(0) uj RhdO , 



(12) 



where a^a is the effective Ha recombination coefficient. 

Taking into account the angular dependence of the density 
profile we can integrate (12) to obtain the Ha emission of the 
cometary heads : 
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nR 



2 a Ha 
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1 + 



4F* 
6 



3/2 



(13) 

The term in curly brackets is the fraction of incident ionizing 
photons that are absorbed in the photoevaporating flow before 
reaching the ionization front. It is only this fraction of the in- 
cident photons that are reprocessed into Ha radiation. If we 
know the size of the knot R/, and the Ha emission Sua we can 
use equation (13) to estimate the stellar flux. 

2.2. A model for the cometary tail 

For the cometary tail, we consider the problem of a cylinder 
of neutral material behind the cometary head being photoion- 
ized by the diffuse flux of the nebula. A schematic diagram of 
this configuration is shown in Figure ||. If we assume that the 
radius of the cometary head is much smaller than the distance 
to the source we can consider the shadow region to have po- 
lar symmetry. This cylindrical shadow does not receive direct 
stellar radiation. We therefore have : 



F d = F, + uj R, n, a B , 



(14) 



where F d is the diffuse ionizing flux of the surroundings, F t is 
the incident flux at the tail surface and the second term on the 
right represents the absorptions in the photoevaporated wind of 
the cometary tail. 

The incident flux at the neutral surface of the tail is related to 
the density n, at the base of the cylindrical wind through : 



F t = n, a . 



Defining the parameter : 



6: 



c? 



uj a B Rt 



and subsituting equations (15-16) in (14) one obtains : 



Fj 



(15) 



(16) 



(17) 



In order to obtain the total number of Ha photons emitted by 
the tails of the cometary knots, we assume that the Ha emitting 
region has a cylindrical section with a thickness h = u>R,. The 
number of Ha photons emitted per unit time and per unit length 
of this cylinder is : 

-^j = 2 ir R, uj R, n] au a , (18) 
and as a function of the parameter £ f the Ha emission is : 

% = 2^,^f, (19) 
Al a B £, 

where Al is a unit of length along the cylindrical tail. If we 
combine equations (17) and (19) we can calculate the diffuse 
ionizing flux as a function of the Ha emission per unit length 
and the radius of the cylinder : 



5h q /A/ a B 
td = ~S — d + 

2 7T R, ana 



Sho/AI a B 
2 it R, a Ha 



1/2 



(20) 



2.3. Diffuse ionizing field inside an H II region 



As well as the radial ionizing radiation field from the Helix 
central star, there will also be a diffuse ionizing radiation field, 
principally due to ground level recombinations of hydrogen in 
the nebula, plus smaller contributions from helium recombina- 
tions and the scattering of stellar radiation by dust grains. For 
the purposes of calculating the properties of the knot tails, the 
important quantity is the lateral flux of the diffuse field, which 
is the flux incident on the surface of an opaque, radially aligned, 
thin cylinder. 

In the "on-the spot" (OTS) approximation, in which all dif- 
fuse ionizing photons are assumed to be reabsorbed by neutral 
H very close to their point of emission, the diffuse flux, Fj, 
across any opaque surface is related to the stellar flux, F*, at the 
same position by (Henney 2000) 



OTS 



Fd 
F* 



4a B n 



(21) 



where ct\ and a B are, respectively, the H recombination rates 
to the ground level and to all excited levels and k = adjo*, 
where dd is the mean photoionization cross-section averaged 
over the diffuse ionizing spectrum and ct* is the same quantity 
averaged over the stellar spectrum. Henney (2000) presents de- 
tailed calculations, in which the OTS assumption is relaxed, for 
the lateral diffuse flux in two simplified geometries: a classical 
filled-sphere homogeneous Stromgren H II region, and a hol- 
low cavity H II region with the ionized gas concentrated in a 
thin spherical shell. It is found that in both cases the value of 
0(= -p-) is much smaller than /3 ts at small radii. It is because 
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F d 




FIG. 2. — Schematic diagram of the cometary tail. F t is the incident flux at the surface of the neutral cylinder, Fj is the diffuse flux (produced by the surrounding 
nebula), and R t and n t are the radius and the density of the neutral tail (respectively). 



at small radi, gets huge (r~ 2 ) but Fj, which is proportional to 
the recombination rate times a path length, plateaus to a more 
or less constant value. In the thin-shell case, (3 remains at a low 
value throughout the interior of the cavity, only becoming com- 
parable to /3ots a t the position of the shell. In the filled-sphere 
case, on the other hand, for the relatively high values of K ap- 
propriate for the Helix knots (see below), j3 rises to ~ /3ots at a 
fractional radius of ~ 0.5 and is roughly constant thereafter. 

The stellar ionizing radiation field is much harder than the 
diffuse field, while the photoionization cross-section declines 
rapidly with frequency above the ionization threshold. Hence, 
the ratio of mean cross-sections, k, for the diffuse and stellar 
fields is larger than unity. Including the H and He recombi- 
nation spectrum gives a a ~ 0.75cto, where <tq is the threshold 
cross-section and He/H = 0.13 (O'Dell 1998; Henry, Kwitter 
& Dufour 1999) has been assumed. 

The temperature of the central star of the Helix has been mea- 
sured by the H/3 Zanstra method to be ~ 1.2 x 10 5 K (Gorny, 
Stasihska & Tylenda 1997). A further constraint on the stellar 
spectrum is the observation (O'Dell 1998) that the He" 1 " 1 " zone 
in the nebula is roughly half the radius of the He + zone, imply- 
ing that roughly 10% of the ionizing photons have frequencies 
higher than the He + ionization limit. Model atmospheres for 
compact hot stars (Rauch 1997) that are consistent with these 
constraints have photon spectra L v jhv that are flat or rising 
between the H and He + ionization limits, implying a value of 
<7* ~ 0.17oo. 

Hence, k ~ 4.5 is appropriate for the Helix knots, which 
implies /3ots — 0.033 for an assumed electron temperature 
of \Q A K. The true situation in the Helix nebula is proba- 
bly intermediate between the filled-sphere and the hollow-shell 
case. Although O'Dell (1998) finds that the electron density is 
roughly constant with radius throughout the nebula, the central 
high-ionization "hole" has a higher temperature and hence a 
lower emissivity of diffuse ionizing photons. Furthermore, the 
geometry is more disk-like than spherical, which would tend to 
reduce the intensity of the diffuse field. 

The radial dependence of j3 in the two limiting cases is shown 
in Figure g. This is compared with the observational data in 
section 3.3. 



3. PROPERTIES OF THE COMETARY KNOTS 

We have used the HST images from the data archive at the 
Space Telescope Science Institute to obtain the Ha emission 
from the cometary knots in the Helix Nebula. These images 
are flux-calibrated using the coefficients given in O'Dell & Doi 
1999. The f656n and f658n filters was selected to isolated the 
Ha and [Nil] emission at 658.4 nm. The [Nil] emission is not 
much stronger than Ha, and thus it does not produce an impor- 
tant contamination of the Ha observations (less than « 15%). 

Several methods have been used in order to obtain the dis- 
tance to this PN (Cahn & Kaler 1971; Daub 1982; Cahn, Kaler 
& Stanghellini 1992; Harris et al 1997), leading to distances 
ranging from 120 to 400 pc. We use the value of 213 pc deter- 
mined by Harris et al. (1997) from trigonometric parallax. 

The comparison between the models and the observations is 
done for 26 knots, chosen for being relatively well isolated, and 
covering a range of distances to the central star. 

3.1. Sizes of the knots 

In order to determine the sizes of the knots, we carry out 
aperture photometry with a series of circular diaphragms (cen- 
tred on the knots) of increasing angular radius. If one plots the 
flux within the aperture versus radius, for large enough radii the 
flux has to increase as the square of the radius (due to the pres- 
ence of the bright, surrounding nebular environment). From a 
logarithmic plot it is then straightforward to determine the ra- 
dius R ext = Rh + ujRh (see equations 4-5 and Figure |l]) at which 
the quadratic flux vs. radius dependence first appears. In this 
way one then determines the value of the radius R/, of the neu- 
tral clump. 

Applying this method to the 26 chosen knots we find that 
there is no correlation between the size of the knots and their 
distance to the central star (see Figure Q). However, this could 
be a result of the fact that we only have a few knots, and that 
these knots are not necessarily representative of the true dis- 
tribution of knot sizes (as they were chosen so as to be well 
isolated, see above). We calculate a mean radius for the knots 
of<R h >= (0.68 ±0.07) ". 
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FIG. 3. — Comparison of the Henney (2000) models to the observed values of the diffuse to direct flux ratio (crosses) as a function of the distance from the 
ionization front. The upper horizontal line is the OTS value and the solid line and dashed line give the variation of /3 across the nebula for the filled-sphere and 
hollow-shell case, respectively. 



3.2. Ha intensities of the cometary heads 

Subtracting the background from the fluxes determined with 
the aperture photometry, we obtain the Ha fluxes emitted by the 
knots. One can then use equation (13) to calculate the stellar 
ionizing flux at the position of the successive knots as a func- 
tion of distance from the central star. 

Column 6 of Table 1 gives the resulting Ha photon pro- 
duction rates Sua of the knots, using values for uj = 0.1 and 
c; = 10 km s" 1 (therefore £/, ~ 1.79 x 10 10 cm" 2 s" 1 for a knot 
with the mean radius < Ri, >= 0.68 ")■ Figure [| shows Sh q vs. 
Di, (where D/, is the projected distance between the source and 
the clumps). 

Figure a also shows the Sua vs. D/, predicted from equa- 
tion (13) for different values of 5* (= 47rD/, 2 /%,). The predicted 
curves should represent upper envelopes of the observed points, 
as the real distances between the clumps and the source are 
larger than the observed, projected distances. From a compar- 
ison of the predicted curves with the values measured for the 
clumps, we see that the observations can best be fitted with a 
stellar ionizing photon rate of 



5* w 5 x 10 45 s" 1 



(22) 



If the knots were distributed isotropically around the central 
star, then the median knot would have a true distance ~ 1.15 
times greater than its projected distance. However, the lower 
envelope of the observed knot brightness distribution suggests 
that all knots lie within ±45° of the plane of the sky, consistent 
with a ring-like spatial distribution. 

According to Osterbrock (1989) 5* can be derived from the 
H/3 luminosity : 



a eS 



(23) 



If we use the value of the H/3 luminosity of the Helix Nebula 
obtained by O'Dell (1998) we have a value for the stellar ion- 
izing photon rate of 5.25 x 10 45 s" 1 . As we can see the match 



with our fit is excellent. 

3.3. Ha intensities of the tails 

In order to obtain the Ha emission from the tails we integrate 
the observed emission over a rectangular area covering the re- 
gion in which the tails are clearly detected. We also integrate 
the emission in two small adjacent rectangular areas in order 
to determine the intensity of the nebular background, which we 
subtract from the emission of the box containing the cometary 
tail. 

Knowing the Ha emission we can then calculate the diffuse 
flux with equation (20) for the 26 chosen knots. In this way we 
calculate the diffuse flux for different distances to the central 
star (see column 7 of Table 1). As we have also computed the 
direct stellar flux from the emission of the heads of the knots 
(section 3.2), we can then calculate the diffuse-to-direct stellar 
flux ratios, and compare these values with the ones predicted 
from the (Henney 2000) models described in section 2.3. This 
comparison is shown in Figure [| where we see that the ob- 
served points fall into the region between the curves delimited 
by the two models, indicating that the case of the Helix nebula 
is intermediate between the homogeneus sphere and the "thin 
shell" cases. 

3.4. Knot masses and evaporation rates 

From the analysis it is also possible to derive knot masses 
and evaporation rates. Under the assumption that the knots are 
accelerating and have an exponential density profile, one can 
derive their mass by 



M- 



1-1 

7T 



16F c (0)mci 



R 3 C , 



(24) 



see Mellema et al. (1998), equation (32). In this equation m 
is the average mass per atom or ion, for wich we take 1.3toh, 
R c is the radius of the clump, and c n is the isothermal velocity 
of sound in the neutral gas for which we take 1 km s" 1 (corre- 
sponding to 150 K). Applying this to the measured knots gives 
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Table 1 

Location and fluxes of the cometary knots 
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FIG. 4. — The angular radii of the knots as a function of projected angular distance to the central star. The crosses indicate the angular sizes determined for 
individual knots. The horizontal, continuous line indicates the mean value of knot radius, and the dashed horizontal lines give the value of the dispersion from the 
mean value. The oblique dotted line gives the best fit to the observed data. 
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FIG. 5. — Model fits to the Ha photon rates observed for the knots (as a function of distance from the central star) obtained for different stellar ionizing photon 
rates. The predicted curves should represent upper envelopes of the observed points. 



an average mass of 9 x 10" 6 M© with knot masses ranging from 
6 x 10" 6 M Q to 1 .8 x 10" 5 M©. In the process one has to calcu- 
late the fraction of the ionizing flux which actually reaches the 
heads of the knots, F c (0)/F*, which comes out to be about 70 
to 80%. 

The mass range found is very close to the one derived by 
O'Dell & Handron (1996) and Meaburn et al. (1998) from dust 
extinction measurements and by Huggins et al. (1992) from CO 
measurements, showing that the sizes and shapes of the knots 
are entirely consistent with the photoevaporation model. 

The model also allows an estimate of the current mass loss 
rate from the heads of the knots, which is given by 



AM 
~d7 



-F c (0)rmrR 2 c 



(25) 



see Mellema et al. (1998), equation (36). This gives an aver- 
age of -2.2 x 10~ 9 Mq yr" 1 (ranging between -4.1 and -1.6 x 
10" 9 M yr" 1 ). 

4. DISCUSSION 

The results of the previous sections show that our photoevap- 
oration model is entirely successful in explaining the current 
situation of the cometary knots in the Helix nebula. One could 
say that we are using the cometary knots as 'probes' for the di- 
rect and diffuse UV radiation in the nebula, under the assump- 
tion that both the heads and tails are photoevaporating. The 
results from the probes are consistent with the results found by 
other means. 

The question which remains unanswered in this approach is 
the origin and final fate of the knots. Starting with the their 
final fate, the photoevaporation models give an estimate of the 
evaporation time. The evaporation time in the "Flux Dominated 



Regime" is f evap « R c c\jc\, so that once a clump start to shrink 
by evaporation, the timescale gets shorter and shorter. The re- 
sults from Mellema et al. (1998) show that the evaporation time 
of a knot is given by 



'evap 



'5 lnd + 7?) 

:/ *' 6 + —&T 



with 



and 



7T V 7T 



2 a B F„R c 



(26) 



(27) 



(28) 



which applied to the measured knots gives an evaporation time 
of 1.1 to 1.2 x 10 4 years. The time for the clumps to move 
from close to the star to their present position is about 6 x 10 3 
years. In the "Flux Dominated Regime" the evaporation time 
is fevap ~ RcC\cl, so that once a clump start to shrink by evap- 
oration, the timescale gets shorter and shorter. There is an un- 
certainty in this number since the sound speed in the neutral 
gas is not well defined. Probably, the knot has a range of tem- 
peratures, depending for instance on how far the molecule dis- 
sociating photons penetrate. The evaporation time is inversely 
proportional to the temperature, the values quoted being valid 
for 150 K. Huggins et al. (1992) cite a temperature of 25 K for 
the CO gas in the knots. 

One consequence of the photoevaporation is that gas is fed 
into the region surrounding the knots. Interestingly enough the 
Helix Nebula seems to be special among PNe in that its 'central 
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cavity' is not empty but is at least partly filled with high ioniza- 
tion gas producing detectable amounts of [OIII] lines, see for 
example O'Dell (1998) and Henry et al. (1999). The density 
of this gas is derived to be around 50 cm" 3 , perhaps less. We 
conjecture that this material was injected into the cavity by the 
photoevaporating clouds. Taking a radius of 200" for the cav- 
ity, and assuming a spherical shape (probably an overestimate 
of the volume and hence of the mass), one finds that the evap- 
oration of 3000 to 5000 knots of about 10" 5 M Q can supply 
this amount of gas. This equals the mass and current number 
of knots in the Helix. Possibly the gas filling the central cavity 
is from already evaporated knots, or from the current knots. In 
Mellema et al. (1998) it was shown that 50% of the mass of a 
photoevaporating clump is lost during the first 'collapse phase'. 
Since the current knots apparently are in the 'cometary phase' 
they might easily have already lost half of their initial mass. 

The question of the origin of the knots can be divided into 
two questions: "Why do they have a cometary shape?" and 
"What physical process is responsible for their existence?". 
Three models exist to explain their shapes. Dyson et al. (1993) 
showed that collisions between a supersonic wind and a clump 
produce short stumpy tails, and that only the interaction be- 
tween a subsonic wind and a subsonic flow from a clump can 
produce long, comet-like tails. When applying this theorem to 
our model we run into the problem that the shapes derived in 
Dyson et al. (1993) are the shapes of the contact discontinu- 
ity between the flow from the clump and the wind, which we 
do not trace in our photoionization description. The dynamical 
models in Mellema et al. (1998) contain both the photoioniza- 
tion processes and the interaction of the photoevaporation flow 
with the environment. In those models the shape of the contact 
discontinuity is far from comet-like, as predicted by Dyson et 
al. (1993); but at the same time the contact discontinuity is not a 
region which produces a lot of emission, as the ionization front 
outshines it by many factors. However, in that model the envi- 
ronment does not exert any ram pressure on the clump and its 
photoevaporation flow. If the ram pressure in the environment is 
high enough it will overwhelm the photoevaporation flow, and 
dominate the dynamics. Observations show that the environ- 
ment of the cometary knots has a density of about 50 cm" 3 and 
a temperature of perhaps 20,000 K. If the environment is flow- 
ing subsonically, this implies pressures of order 10~ 10 dyn cm" 2 . 
The photoevaporation pressure is given by 2F c mc; (see Mellema 
et al. 1998), which gives values of about 10" 8 dyn cm" 2 . Clearly 
the environment will not be able to overwhelm the photoevap- 
oration flow from the knots. 

We cannot rule out the possibility that some time in the 
past the cometary knots were shaped into their cometary shape 
through a wind-clump interaction. In fact, below we will argue 
for such a scenario to explain the origin of the knots. 

Canto et al. (1998) suggested that long tails can be formed 
behind clumps opaque to ionizing UV photons. The region be- 
hind the clump does not see any direct ionizing photons from 
the star, just the diffuse UV radiation field, which in general 
leads to denser, partly neutral tails forming behind the clumps. 
However, these tails will be initially the same density as the 
environment, fill up and recombine, striving towards pressure 
equilibrium with the ionized environment. The tails of the He- 
lix knots are overpressured compared to their environment, a 
situation which does not occur in the scenario put forward in 
Canto et al. (1998). The overpressured tails suggest that they 
were formed under different circumstances from the ones in 
which they are now. 



Burkert & O'Dell (1998) proposed that the knots are shaped 
by Lya photon radiation pressure. This was partly motivated 
by the fact that the evaporation flows appear to have exponen- 
tial brightness profiles. These autors found that this mechanism 
does not work in the cometary knots of the Helix nebula. A 
similar model was proposed for the Proplyds in the Orion neb- 
ula (O'Dell 1998). Henney & Arthur (1998) showed that the 
exponential brightness profiles are not inconsistent with photo- 
evaporation flows and that the radiation pressure from Lya is 
at least one order of magnitude too low to be significant. Their 
arguments also hold for the Helix Nebula cometary knots and 
we will not repeat them here. The conclusion is that none of the 
three models can explain the origin of the long cometary tails 
under the current circumstances. 

When we consider the origin of the knots there are basically 
two options, the knots are either primordial or the result of in- 
stabilities. Meaburn et al. (1998) and Young et al. (1999) favour 
the primordial model, the main arguments being that the CO gas 
outside of the main nebula has a very clumpy distribution, and 
that the velocity distribution of the cometary knots follows that 
of the main CO emission, albeit at a lower expansion velocity. 
Within this model it is possible that the tails were shaped as the 
knots were run over by the main part of nebula. During this 
phase they would be surrounded and eroded by subsonically 
flowing gas, which according to the Dyson et al. (1993) model 
would produce long thin tails. Also, the formation of the tails 
as ionization shadows behind the clumps as proposed by Canto 
et al. (1998) would in this case produce long, dense tails. After 
the main nebular shell has passed by, the knots find themselves 
fully exposed to the ionizing stellar radiation and the tails to 
the diffuse UV radiation field and both start photoevaporating, 
which is the stage in which we see them now. 

O'Dell & Handron (1996) prefer the model in which the 
knots are formed as Rayleigh-Taylor instabilities on the inside 
of the swept up main nebula, proposed originally by Capri- 
otti (1973). The numerical model of R-T instabilities presented 
in figure 4 in O'Dell & Burkert (1997) seems to support this 
model; it shows elongated structures pointing radially to the 
centre of the nebula. However, the much higher resolution mod- 
els of Walder & Folini (1998) show the R-T instabilities to be 
much more chaotic. The fingers in their model need a lot of 
processing before they look anything like the cometary knots, 
and one would also expect much more variation in knot prop- 
erties if they were formed in such instabilities. In all we prefer 
the first option of primoridal knots. 

5. CONCLUSIONS 

In this paper we propose an analytical model for the heads of 
the cometary globules in the Helix Nebula, in which the emis- 
sion is assumed to come from a flow which is being photoe- 
vaporated from the surface of a neutral clump. Using the Ha 
emission obtained from HST images of the Helix Nebula (ex- 
tracted from the HST archive, see also O'Dell & Handron 1996; 
O'Dell & Burkert 1997) and the predictions from our analytic 
model, we can calculate the ionizing stellar photon flux at the 
positions of the knots. 

We find that the knot brightnesses are fully consistent with 
the ionizing photon luminosity of S+ = 5.25 x 10 45 s" 1 for the 
central star, deduced from the total Hf3 flux of the nebula. This 
value, however, is somewhat uncertain as in the present work 
we have made no attempt to substract the [N II] 6583 flux from 
the f656n filter. Given the rather strong [N II] emission of some 
of the Helix Knots, our Ha flux determinations could in some 
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cases lie up to «30 The angular size of the knots has a value of 
< R h >= (0.68 ± 0.07) " (where R h is the radius of the neutral 
clump forming the knot, see section 3.1). From our small sam- 
ple of 26 knots, we find no correlation between R/ t and distance 
from the central star (see Figure Q). 

The photoevaporation model also predicts masses for the 
knots, which come out fully consistent with observationally de- 
rived masses. Their evaporation time is between several thou- 
sands to 10 4 years, dependent on the temperature of the neutral 
gas in the knots. 

In this paper we also present a model for the cometary tails 
as a flow which is being photoevaporated from a neutral, cylin- 
drical "shadow region" behind the neutral clumps which form 
the knots. This flow is photoevaporated by the impinging dif- 
fuse ionizing field, formed in the surrounding nebula. As we 
have done for the heads of the cometary globules, we compare 
the emission of the cometary tails with the predictions from our 
model in order to quantitatively determine the value of the dif- 
fuse ionizing photon flux at the positions of the globules. 



Knowing the direct stellar radiation and the diffuse flux at 
different distances from the central star (i.e., at the positions 
of the cometary globules), we can then calculate the diffuse-to- 
direct ionizing flux ratio (Fd/F*) as a function of distance from 
the central star. We compare the data with simple models for 
the diffuse ionizing field (Henney 2000) and conclude that the 
diffuse field is intermediate in magnitude between that expected 
from a filled and shell-like geometry, which may be due to the 
disk-like nature of the Helix nebula. 
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